Inclusive B decays from resummed perturbation theory 
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I review the recent progress in computing inclusive B decay widths and spectra and its implica- 
tions on the interpretation of measurements from the B factories. I discuss the inclusive charmless 
semileptonic decay, B —* X U W, which provides the most robust determination of the CKM param- 
eter |V u b|, and the rare radiative decay, B — > X s "/, which constrains flavor violation beyond the 
Standard Model. I demonstrate that precise predictions for the experimentally-relevant branching 
fractions can be derived from resummed perturbation theory and explain the way in which the 
rcsummation further provides guidance in parametrizing non-perturbative Fermi-motion effects. 
Finally I address the comparison between theory and data and discuss future prospects. 
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I. INTRODUCTION 

Inclusive B decay measurements have a central role 
in flavour physics. Inclusive semileptonic decays [l| 
B — > X CA1 lv provide the most accurate determinations of 
the CKM parameters |V^b| an d IKihL an d thus strongly 
constrain the Unitarity Triangle [2, 3]. This constraint is 
particularly important: being based on tree-level Weak 
decays it is totally insensitive to physics beyond the Stan- 
dard Model (SM). This should be contrasted for example 
with the measurement of the CP violating phase sin(2/3), 
which is sensitive to loop effects through B° — B° mix- 
ing, and can therefore be influenced by physics beyond 
the SM. The potential to discover new physics crucially 
depends on our ability to make precise comparison be- 
tween the two @, i, H i, H . 

Another major avenue in flavor physics is the use of 
rare decays to constrain flavor-changing neutral currents 
beyond the SM. Here the most stringent constraint is 
based on the inclusive Branching Fraction (BF) of radia- 
tive B decays, B — > X s j, see [3]. 

Our ability to exploit the potential the B factories de- 
pends on the accuracy at which we can compute the 
corresponding BF's. Since only part of the phase space 
is experimentally accessible, theoretical calculations are 
needed not only for the total widths, but also for the 
spectra. 

In B — > X s j and in B — > X u lv most events are charac- 
terized by jet-like kinematics, where the invariant mass 
of the hadronic system mx is small, much smaller than 
the energy 0(mb) that is released in the decay. More- 
over, the measurements are restricted to this kinematic 
region. Therefore, theoretical understanding of the limit 
rax -C rrib becomes the key to the interpretation of the 
measurements. 

A major part of this talk is devoted to the re- 
cent progress in our ability to compute inclusive de- 
cay spectra using resummed QCD perturbation the- 
ory @,H,[l(i[n]>[3- The most precise predictions for the 
experimentally-relevant BF's in B — > X s j [l(J E3] and 
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FIG. 1: HFAG world average |V u b| using DGE: summary of all 
available inclusive measurements of B — * X U W with different 
kinematic cuts on Mx, q 2 , Ei etc. In each case the partial BF 
was translated to a value of |Vub| using DGE. 



in B — > X u lv [TlJ (see Fig. [T]) are now based on calcu- 
lations using resummed perturbation theory — not any- 
more on a parametrization of experimental data using a 
"shape function" . This breakthrough was achieved using 
Dressed Gluon Exponentiation (DGE) [13j] . a resumma- 
tion method that fully uses the inherent infrared-safety 
of the on-shell decay spectrum. I will explain the con- 
ceptual differences between DGE and the "shape func- 
tion" approach, show that definite predictions can be de- 



rived from perturbation theory despite its divergent na- 
ture, and finally explain the way in which the resumma- 
tion provides guidance in parametrizing non-perturbative 
Fermi-motion effects. I will also discuss the determina- 
tion of the B — > X s -f BF and conclude by addressing the 
present theoretical uncertainty, the comparison between 
theory and data and future prospects. 



II. BREAKTHROUGH IN COMPUTING 
DECAY SPECTRA 

The key to computing inclusive decay spectra is the 
understanding of the Sudakov region, where the hadronic 
system in the final state is jet-like, namely it is character- 
ized by a large hierarchy between the two lightcone mo- 
mentum components, P + <C P~ , where P ± = ExT\Px\ 
and P~ is typically of the order of the heavy-quark mass. 

It has been understood long ago that the leading con- 
tribution to the spectrum for small P + involves the mo- 
mentum distribution of the b quark in the meson [T^.[l5|. 
Indeed, if we imagine a hypothetical situation where the 
b quark (of mass mj) is on-shell inside the meson — so it 
does not interact with the light-degrees-of-freedom car- 
rying the residual energy A = Mb — nib — we would 
find that the spectrum, at any order, has support only 
for p + = P + — A > 0. In contrast, the physical spectrum 
has support for P + > 0. Fig. [2] demonstrates these dif- 
ferent support properties in the case of B — > X s j where 
P~ = M B and P+ = M B - 2E 7 . It is therefore clear 
that the spectrum 1 at small P + is dictated by the Fermi 
motion of the b quark inside the meson. 

The conventional approach to describing the spectrum 
has been based on the analogy with deep inelastic struc- 
ture functions: the non-perturbative momentum distri- 
bution function is parametrized and is convoluted with 
a hard coefficient function that is computed with an in- 
frared cutoff. In this approach the ansatz for the momen- 
tum distribution function essentially dictates the shape 
of the spectrum in the entire peak region, hence the name 
"shape function" . This is unfortunate, as there is no the- 
oretical guidance whatsoever as to how to parametrize 
this function. All that can be currently deduced from 
fits to data are the center of the "shape function" (the b 
quark mass in some definition) and its width. Even with 
much more precise data in the peak region, the tail of the 
"shape function" would not be well constrained. Thus, 
the prospects of improving the precision in this approach 
(beyond what has already been done [HI, are very 
limited. 



1 The physical spectrum near the exclusive limit depends also on 
the hadronic structure of the final-state jet (e.g. the B — > Xs'y 
spectrum contains the B — * K*^ peak) but inclusive obscrv- 
ables such as the partial branching fractions and the moments 
(defined over a sufficiently large phase space) are insensitive to 
these details. 
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FIG. 2: The photon-energy spectrum in B — > X s j (for 
the electromagnetic dipole operators O7) as computed by re- 
summed perturbation theory (DGE) compared with fixed- 
order perturbation theory at NLO, NNLO [IH, Q3] and higher 
order corrections in the large-/3o limit . The endpoint of the 
on-shell decay spectrum is /2 (arrow) where the quark 
pole mass is v = 4.87 GeV. The endpoint of the physical 
spectrum is Mb /2 (end of the E-y axis) where the meson mass 
is M B = 5.28 GeV. 



A more promising avenue is that of Dressed Gluon Ex- 
ponentiation (DGE) [l3l ] where the first approximation 
to the spectrum is based on the resummed calculation 
of the on-shell decay spectrum; non-perturbative effects 
enter only as power corrections M. DGE is a general 
resummation formalism [3, [2(J, H3] for inclusive dis- 
tributions near a kinematic threshold 2 . It goes beyond 
the standard Sudakov resummation framework by incor- 
porating renormalon resummation in the calculation of 
the exponent. This has proven effective [2(J [HJ in ex- 
tending the range of applicability of perturbation theory 
nearer to threshold and in identifying the relevant non- 
perturbative corrections. 

The impact of resummation on the calculation of the 
B — ► X s j spectrum is demonstrated in Fig. O In fixed- 
order perturbation theory one finds very large correc- 
tions. The parametrically-leading corrections near the 
partonic endpoint E 7 — > rrib/2 (p + — > 0) are Sudakov 
logarithms. However, conventional Sudakov resumma- 
tion with a fixed logarithmic accuracy yields a divergent 
expansion (see Fig. 1 in 23]) as it misses out important 
running-coupling corrections. In contrast the DGE re- 
summed spectrum yields a stable result, see Fig. [5] 

A remarkable difference between the fixed-order re- 



2 Here the threshold corresponds to the phase space limit where 
the small lightcone momentum component vanishes, p+ — » 0. 
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suit and the resummed one is in the support proper- 
ties. While the perturbative endpoint, at any order, is 
E 1 = rrib/2 (i.e. p + = 0), the resummed result extends 
into the non-perturbative regime and tends to zero near 
the physical endpoint E 1 = M B /2 (i.e. P+ = 0). Thus, 
resummation makes a qualitative difference. In the fol- 
lowing we briefly describe what is done to arrive at this 



result, and then proceed to explain how Fermi-motion 
effects arc taken into account as power corrections. 

When applying perturbation theory to inclusive decay 
spectra it proves useful to consider the moments with 
respect to the ratio between the lightcone momentum 
components, e.g. in the semileptonic decay we define: 
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dT b ^ x - w (p+,p-, Ei 
dp + dp~ dEi 



(1) 



where the partonic lightcone momentum components p 
are related to the hadronic ones by: p — P — A, where 
A = Mb — rrib is the energy of the light-degrees-of- 
freedom in the meson. Note that we consider the mo- 
ments of the fully differential width [Tlj : e.g. in the 
semileptonic case the moments remain differential with 
respect to the large lightcone component p~ as well as 
the lepton energy E\. This is essential for performing 
soft gluon resummation. 

Owing to their inclusive nature the moments |T]) are 
infrared safe: for any moment N there is an exact can- 
cellation of infrared singularities at p + — » between real 
and virtual diagrams. Thus, instead of considering di- 
rectly the singular limit p + — > 0, one considers the large- 
TV limit. In this limit one identifies three characteristic 
scales, and the moments factorize [E [IE S3, SH to all 
orders as follows: 

dp- dE t 1 ' 

H{p-, E^ J(p-/VN, n) S b (p-/N, ft) +0(1/N), 

" v ' 

Sud(p-, N) 

where the factorization-scale (/i) dependence cancels ex- 
actly in the product in the Sudakov factor Sud(p _ , N). 
Factorization facilitates the resummation of Sudakov log- 
arithms @, EE Hi HE SUM Hi (see also [MSESEHJ), 



the corrections that dominate the dynamics at large N. 
These originate in two distinct regions of phase space: the 
final-state jet J(p~ /y/~N, /i) and the initial-state quark- 
distribution function Sb (p~ /N, fx) . The latter subprocess 
depends on the softest scale, p~ /N, and is therefore the 
first place where non-perturbative corrections should be 
included. Indeed, the meson decay is described by 

arfr x ^(p-, go = m 

dp- dEi 1 1 

H(p-, E^ J{p~/s/N, n) S B (p-/N, n) + 0(1/N), 

where Sg(p-/N,/i) stands for the quark distribution in 
the meson. Owing to the Fermi motion, this function 
differs from its perturbative counterpart, the quark dis- 
tribution in an on-shell quark Sb(p~ /N, /i) , by power cor- 
rections (AN/p~) k , which can be resummed into a new 
non-perturbative function, 

S B (p-/N,») = S b (p-/N,n) x F(p-/N). (4) 

Note that T(p~ /N), in contrast to S§(p- /N, (i) , is in- 
dependent of the factorization scale /x. In Eq. (J5J) below 
we shall parametrize T{p~ /N) based on the infrared sen- 
sitivity exposed by the resummation. Finally, the spec- 
trum is given by an inverse Mellin transform to ([1]) : 
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-N 
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(p-,Et 



dp dEi 



(5) 



p±=P±-A 



where the integration contour runs parallel to the imag- 
inary axis, to the right of the singularities of the inte- 
grand. 

We are now in a position to understand some of the 
fundamental differences between the approach described 



above and the one based on a "shape function" : 

• When considering the on-shell decay in moment 
space for N ^> 1 we begin by assuming that 
p~ (i.e. mb) is sufficiently large so that even 
the soft scale characterizing the quark distribu- 
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tion function, p~ /N, is in the perturbative regime: 
p~ /N ^> A. In this limit the on-shell decay spec- 
trum is computable and provides a systematic ap- 
proximation to the meson decay spectrum. This 
should be contrasted with the "shape function" 
approach where one works in momentum space 
and considers directly the situation where p + = 
0(A), so the "shape function" is treated as non- 
perturbative from the beginning. 

• Consequently, we can compute more and 
parametrize less. We compute the quark dis- 
tribution in an on-shell quark, Sb(p~/N, /i) , which 
accounts for the radiation off the heavy quark. 
This radiation puts the heavy quark slightly off its 
mass shell. Thus, the quark decays off its mass 
shell, despite the fact that it was initially assumed 
on-shell. Of course, physically there is no unique 



distinction between this acquired virtuality and 
the "primordial" virtuality, which is "a property of 
the bound state" . This corresponds to the renor- 
malon ambiguity in the calculation of the quark 
distribution function, which serves as a probe of 
the non-perturbative dynamics. Thus, eventually 
we need to parametrize the power corrections 
T(jT /N) that make for the difference between the 
quark distribution in the meson S§(p~/N,ii) and 
that in an on-shell quark. In contrast, the "shape 
function" approach is based on parametrizing the 
entire function Sg(p~~ /N, fu . 

The functions j(p~/y/N, pi) and Shfo~JN^j) in © 
satisfy Sudakov evolution equations [ij l27l. l28l. 123. l30l. I3U 
[HJ whose general, all-order solution can be formulated 
as a scheme-invariant Borel sum [IE [II 11,11: 



Sud(p , N) = exp 



fCV f°° du /A_ 
I A) Jo u \p~ 



2u r 



B s (u)T(-2u) 



T(N) 



1 



T(N-2u) 
Bj(u)T(-u) 



r(i - 2u) 

T(N) 



r(i 



(6) 



where B$ (u) and B j (u) are the Borel representations of 
the Sudakov anomalous dimensions of the quark distri- 
bution and the jet function, respectively. The Sudakov 
exponent has renormalon singularities at integer and half 
integer values of u, except where Bg,j{u) vanish. The 
corresponding ambiguities, whose magnitude is deter- 
mined by the residues of the poles in ([6]), are enhanced 
at large N by powers of N. They indicate the presence 
of non-perturbative power corrections with a similar iV 
dependence. These power corrections exponentiate to- 
gether with the logarithms. 

The anomalous dimension functions Bs{u) and Bj{u) 
are both known [9] to NNLO, 0(u 2 ). This facilitates Su- 
dakov resummation with next-to-next-to-leading loga- 
rithmic accuracy (see also P, EM)- In the DGE 
approach one takes a crucial step further [2(| Hl|: in- 
stead of expanding ([6]) and truncating the expansion at 
a given logarithmic accuracy, the Borel integral in ^ - 
the "renormalon sum" — is performed in the Principal 
Value (PV) prescription, incorporating information on 
the behavior of the integrand away from the origin [l(| • 
This information includes, in particular, the pattern of 
zeros of Bg,j{u), which one identifies using the large- 
(3q limit: a zero implies the absence of a corresponding 
renormalon ambiguity in the exponent, a property we call 
"infrared safety at the power level" [l3| . 

The renormalon ambiguities probe the infrared sen- 
sitivity and can therefore be used to parametrize non- 
perturbative power corrections. Let us consider the 



power corrections associated with the quark distribution 
function. Upon choosing the PV prescription in the ex- 
ponent in ([6]) we essentially define the quark distribution 
in an on-shell quark at the power level. In Eq. f4|, 
this amounts to a specific separation between perturba- 
tive and non-perturbative contributions to the (unique!) 
quark distribution in the meson. Operationally, promot- 
ing the resummed partonic calculation ([2]) into a pre- 
diction for the meson decay, ^ to ([5]), amounts to the 
following replacement: 

Sud(p-,iV) — » Sud(p-, A0| py x Fpv( P -/N) (7) 

where J r py(p~ /N) is constructed as a sum over the 
residues of the Borel integral ((6]) [l2j : 

^ PV (p-/iV) = exp(^7r/pv X (8) 

where we introduced a single dimensionless coefficient 
/pv controlling the magnitude of the power corrections. 
On general grounds fpy is expected to be of order 1: 
non-perturbative corrections are typically of the order of 
the ambiguity. Note that the sum in ([5]) excludes the two 
leading powers k = 1 and k = 2. The leading renormalon 
(k = 1) corresponds to the ambiguity in defining the 
pole mass and it cancels exactly [8J in the physical spec- 
trum ([5]). The second potential renormalon (k = 2) does 
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FIG. 3: Theoretical prediction for the partial BF as a func- 
tion of the cut on the photon energy, _E 7 > Eo, and the first 
two moments defined over the same range ((5.3) and (5.4) in 
[lOl|'). The moments are compared with available data. No 
fits are done here. The effect of potential power corrections is 
illustrated by the difference between two curves: full line with 
C3/2 = 1; /pv = and dotdashed one: C3/2 = 6.2; fpy = 0.3. 

not appear Q : this is an example of infrared safety at the 
power level. Consequently, the leading power correction 
corresponds to k = 3 in ([8]): it scales as (NA/p~) 3 . This 
high power makes the on-shell decay spectrum a good 
approximation to the physical meson decay. 

It should be emphasized that although formally all 



power corrections in ([5]) contribute at sufficiently large N, 
because of the factorial suppression only the first few 
powers are important. This facilitates the practical use 
of ([5]) despite the fact the behavior of B$(k/2) for k > 3 
is not known 3 . We describe the Fermi- motion effects us- 
ing two parameters [l2T ] in ([8]): the k — 3 renormalon 
residue Bs(u = 3/2) = — 0.23366 C 3 / 2 and the overall 
power correction coefficient fpy. As shown in Fig. [21 
in the absence of power corrections, the support proper- 
ties of the DGE spectrum are close, but not identical, to 
the physical ones. One can therefore constrain the pa- 
rameters (C3/2, /pv) [12] using the support properties. 
Reassuringly, one finds that the typical magnitude of fpy 
is <J 1, as expected on general grounds. 

Fig. [3] presents theoretical predictions for the B — » X s j 
BF as well as the average and the variance of the photon 
energy as a function of a cut E 1 > E . The figure illus- 
trates the potential effect of power corrections. The total 
theoretical uncertainty is shown by the green band. A 
detailed analysis of the individual sources of uncertainty 
reveals that for any cut Eq < 2.2 GeV, the uncertainty 
in the average energy is dominated by the parametric er- 
ror assumed for the b quark mass (m" s = 4.20 ± 0.04); 
for Eo > 2.2 GeV it is dominated by power corrections. 
The uncertainty in the variance is dominated by power 
corrections over the whole range Eo > 1.8 GeV. The con- 
clusion is clear: a fit of the computed spectrum 4 to the 
measured moments can yield an accurate determination 
of the b quark mass and the power corrections. Fig.[3]also 
presents the available data for the average energy and the 
variance. The data agrees well with the calculations over 
the entire range of cuts. This comparison suggests that 
the data is already precise enough to constrain the three 
parameters (to;,, C3/2, fpy)- 

Let us return now to the determination of |V u b| from 
inclusive B — » X U W measurements. The HFAG compi- 
lation in Fig. [1] summarizes the extracted values of |V u b| 
from all available measurements based on different kine- 
matics cuts. The crucial ingredient in this determina- 
tion is the calculation of i? C ut, the partial BF within the 
region of measurement over the total b — > u BF. The 
consistency of |Vub| between the different measurements 
provides a good test of the theory predictions as well as 
the experimental analysis. Nevertheless, some more fo- 
cused analysis can directly test the error estimates and 
shed light on specific issues such as the relevance of Weak 
annihilation. Since the DGE calculation 5 is available at 
the level of the triple differential rate, any conceivable cut 
can be readily implemented. The calculation can thus be 
used to check the stability of the extracted |V u b| as a 
function of the cut and, eventually to optimize the cuts. 



3 Bg(u) is known analytically only in the large— /3n limit. 

4 A C++ code is available at: 

http : //www . hep . phy . cam . ac . uk/~andersen/BDK/ 

5 A C++ code is available at: 

http: //www. hep. phy . cam. ac .uk/~andersen/BDK/B2U/ 
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We note that the present theoretical error in determin- 
ing |Kib| is comparable to the experimental one. There- 
fore, making good use of future data from the B fac- 
tories requires more precise theoretical predictions. Let 
us consider for example the case of the Mx < 1.7 GeV 
cut, where a detailed theoretical-error analysis has been 
done in [ill ]. The three dominant sources of uncer- 
tainty on R cu t(Mx < 1.7) are: parametric error (~ 7%), 
higher-order corrections (~ 6%), and power corrections 
(<; 3%). All three can be reduced in the foreseeable fu- 
ture: the parametric uncertainty is dominated by the 
short-distance b quark mass, the higher-order correc- 
tions are dominated by running-coupling effects which 
have been recently computed [35| , and the leading power 
corrections can be deduced from fits to the measured 
B — ► Xg'y moments. Thus, the prospects for improving 
the determination of |Vub| are high. 

III. THE B -> X sl BRANCHING FRACTION 

The B — > X s j BF provides a crucial constraint on 
flavor-changing neutral currents beyond the SM. Exper- 
imental measurements by the B factories are getting more 
precise Q: 

B (B -> A s7 , £ 7 > 1.6 GeV) = (355 ± 26) • 1CT 6 . (9) 

Consequently, a very significant effort has been invested 
in the recent years in computing the SM BF with NNLO 
accuracy (see e.g. [H, H3] and refs. therein). Despite 
the fact that only partial results are available — and, in 

I 



where Ci (/i) are coefficient functions associated with dif- 
ferent operators in the effective Weak Hamiltonian and 
Gij(Eo^) are the corresponding B — > X s ~f matrix el- 
ements. Here we made the important distinction be- 
tween the short distance mass m^ s and the pole mass 
m b , which enter the expression for the width via the op- 



fact, very important ingredients such as the NNLO ma- 
trix element G27 are missing — first numerical estimates 
for the total BF have recently been published [H, H3] : 

B (B -> A s7 , El > 1.6GeV) = {[^ ^ 23)T0- !$ 10 > 

where the result of fl2l includes only (3oa% terms [38[ at 
NNLO while that of [37| includes also partial results and 
estimates of other NNLO terms. Although the two re- 
sults in (fTTTj) are consistent, the theoretical uncertainty is 
substantially different. To have a more complete picture 
it is useful to compare the above estimates to NLO ones. 
Ours is 

B(B —> A s7 , £ 7 > 1.6GeV) = (313 ± 45) • 10" 6 [12] (11) 

while the well known result by Gambino and Misiak from 
2001 reads: 

B(B -> X s j,E 1 > 1.6GeV)= (360 ± 30) TO" 6 [39]. (12) 

The various analysis differ in many ways, and a detailed 
comparison goes beyond the scope of this talk. Here I 
would like to focus on one crucial aspect concerning the 
evaluation of the total BF, namely the way in which the 
b quark mass is treated. I will then briefly discuss the 
cut dependence. 

The dependence of the BF on the b quark mass is 
important since the mass enters the expression for the 
width at the fifth power: 



(13) 



erator and via the phase-space integration, respectively. 
Since the pole mass suffers from an 0(A) renormalon am- 
biguity it cannot be directly used when evaluating the 
expression for the width at fixed order. A similar prob- 
lem is encountered in other decays, for example in the 
case of the charmless semileptonic width one has: 



T(B -> X sl ,E 7 > E ) = \V th V t l\ 2 (rnf\m b )) Z m 3 b £ C^C^) G tJ (E , fi) 



327r 4 l^tbVy [rn b °(m&) 



/oM + /i(A0— + /2OO 



Mm) 



r b 



X u lv, E~, > Eq 



G 2 f 



Vnb\ 2 m 5 b 

1927T 3 



1 + Sl- 



.(/*) 



S2(M) 



Mm) 



(14) 
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Since this renormalon ambiguity cancels exactly [IE EH, 
[42| if both the pole mass and the on-shell decay width 
(G u in (TTJ| and F in (fT3")0 are systematically considered 
to all orders, it does not necessarily limit the precision in 
evaluating the BF's. Two main strategies have been pro- 
posed: (1) the use of an alternative mass scheme, which is 
renormalon free; (2) Borel summation of the expansions 
for the pole mass and the decay width [11, 43]. The first 
strategy is simple, as it can be applied at fixed order, 
however, it involves an uncontrolled uncertainty owing 
to the choice of the mass scheme. This strategy was used 
in [13, to obtain the results quoted above. The second 
procedure requires renormalon resummation, but it has 



two important advantages: it deals directly with the most 
important corrections, those associated with the running 
coupling, and it facilitates using additional information 
on the large-order asymptotic behavior of the expansion. 
This approach has been used in the recent calculation of 
the total charmless semileptonic width In the case 
of the radiative decay, however, there is an additional 
complication owing to the renormalization of the oper- 
ators and their mixing. This gives an advantage to the 
fixed-order treatment where the scale dependence can be 
used as an uncertainty measure. Instead of using some 
arbitrary mass scheme one can normalize the radiative 
decay width using the semileptonic width as follows [12] : 



T(B — X sl ) = \V th V*\ 2 x ml 



F/Gi 



3/5 



Rcsummcd: prescription independent 



Fixed Order 



(15) 



where each of the factors is both renormalization-group 
invariant and renormalon-free: the product to^Gu^ 5 can 
be taken directly from the semileptonic decay analysis, 
while the function F/G 3 ^ 5 is simply treated at fixed or- 
der. This new strategy was used in [l2| to obtain the 
results quoted in (TiT]) and (JTUJ) above. It turns out that 
in this formulation, and with only floa 2 terms at NNLO, 
the renormalization scale dependence is still over 10%. 

Our main conclusion is that the theoretical uncer- 
tainty on the total BF is not yet significantly smaller 
than at NLO. One reason for that is the large can- 
cellation between different contributions to the width, 
e.g. between the G77 sector, which is known in full to 
NNLO [TEES El El, and the G 27 one which is not. It 
is quite clear that the uncertainty would indeed reduce 
upon completion of the NNLO calculation. Finally, com- 
paring different strategies to deal with the b quark mass 
is absolutely essential for having a reliable error estimate. 

Let us turn now to the dependence of the partial BF 
on the photon-energy cut. As discussed in Sec. [Til the 
DGE approach opens the way for making quantitative 
predictions for the spectrum. As shown in Fig. [3] there 
is good agreement between the calculation and the mea- 
surements of the first and second moments of the photon 
energy over the entire range of cuts, 1.8 < Eq < 2.3. 
Moreover, assuming similar cuts, the relative importance 
of non-perturbtaive corrections in the partial branching 
fraction is lower than in the moments. This gives us con- 
fidence in the predictions for the partial BF. 



IV. CONCLUSIONS 

Over the past few years we made significant progress 
in computing inclusive B-meson decay spectra. This 



progress was achieved by DGE, which replaces the "shape 
function" by resummation of the perturbative expansion 
and parametrization of the dominant non-perturbativc 
power corrections. This framework provides the most ac- 
curate predictions for the experimentally relevant BF's in 
inclusive decays. It has passed two crucial tests in com- 
parison with data: 

• The consistency between the extracted values of 
|Vub| from all available B — » X u lv measurements 
with a variety of cuts, as summarized in Fig. [TJ 

• The agreement with the measurements of the first 
and second moments of the photon energy in B — > 
X s -f for all cuts, as shown in Fig. [3) 

So far the theoretical predictions for the spectra [IE EH, 
EIJ have been based on resummed perturbation theory for 
the on-shell heavy quark decay, while non-perturbative 
Fermi-motion effects have only been considered as part 
of the error analysis. As discussed in Sec. [TTJ a main ad- 
vantage of DGE is the direct link it makes between the 
perturbative and non-perturbative sides of the problem: 
definite predictions for the parametric form of the power 
corrections emerge from the resummation formalism. By 
using experimental data to determine the power correc- 
tions, the full predictive power of DGE can be put to 
use. 

In both B — > A s 7 and the B — > X u lv the dominant 
theoretical uncertainties (aside from the sensitivity to the 
b quark mass) are still associate with perturbative correc- 
tions. In both cases, there has been significant progress 
on NNLO calculations but major challenges on the way 
to this goal still lie ahead. 

In B — > X s 7, despite tremendous progress in NNLO 
calculations made by several groups, the theoretical un- 
certainty in the total BF is still larger than the exper- 
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imental one. A significant reduction of the uncertainty 
is expected upon completion of the NNLO calculation. 
As in other cases, also here: the measurements are (so 
far) consistent with the SM predictions — compare © 
with PP . 

In B — > X U W, the fully differential width has recently 
been computed to all orders in the large-/3o limit (3f| . By 
matching the DGE spectrum to the newly available fioa 2 s 



result we expect to improve the determination of |Vub| 
from all kinematic cuts. Further reduction in the theo- 
retical uncertainty can be achieved by more accurate de- 
terminations of mt and constraints on power corrections. 
These require dedicated comparison between theory and 
data for both the B — > X s j moments and for partial BF's 
and moments of the charmless semileptonic decay itself. 
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